Comment on 

"Critical behavior of a two-species reaction-diffusion problem" 
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In a recent paper, de Freitas et al. [Phys. Rev. E 61, 6330 (2000)] presented simulational results 
for the critical exponents of the two-species reaction-diffusion system A + B — > 2B and B — > A 
in dimension d = 1. In particular, the correlation length exponent was found as v = 2.21(5) in 
contradiction to the exact relation v — 2/d. In this Comment, the symmetry arguments leading 
to exact critical exponents for the universality class of this reaction-diffusion system are concisely 
reconsidered. 

PACS number(s): 64.60.Ak, 64.60. Ht, 05.40. +j 

In a recent paper, de Freitas et al. [Q presented a Monte Carlo study of the two-species reaction-diffusion system 
A + B — > 2B and B — ► A in dimension d = 1. They reported the values /3 = 0.435(10) and v — 2.21(5) for critical 
exponents of the order parameter and the correlation length, respectively. The measurement of the short time scaling 
exponent 6' |2) seems consistent with the scaling laws 6' = —rj/z and 2/3 = v{d + rj) using the exact value z = 2 of 
the dynamical exponent. The critical reaction-diffusion system above || belongs in general to the universality class 
of directed percolation (DP) processes coupled to a secondary conserved density that I call DP-C Q in analogy to 
the Model C. In the same manner as the universal behavior of the critical dynamics of a relaxing non-conserved order 
parameter near equilibrium (Model A) is changed to Model C by the coupling to a conserved density, DP processes 
are changed to DP-C processes, de Freitas et al. assume equal diffusion constants for both species. Therefore a special 
DP-C process called KSS and identified by Kree et al. [0 several years ago describes the special system studied here. 
In the appendix of their paper, Kree et al. show by means of a Ward identity that the correlation length exponent 
obeys the exact relation v = 2/d. This yields v — 2 in one dimension. Thus the value of v reported by de Freitas et 
al. cannot be correct. Likewise their conjectured simple fractions for the critical exponents have to be rejected. 

Because the arguments leading to exact critical exponents of the DP-C processes are more or less implicit in several 
papers H,^), I will reconsider their derivation in this comment, and show that they all have their roots in particular 
symmetry properties. 

The Langevin dynamics of the DP-C class can be described by the dynamic functional 



J 



dt d d x < 8 



d t + A(r - V 2 + fc) + ^(gs- gs) 



d t c--iV 2 {c + <7s) -7(Vc) 



(1) 



Here s and c are the densities of the percolating agent and the conserved field, respectively. In the case of the 
reaction-diffusion system above, s oc ng and c oc ha + tib, where tia and tib denote the densities of the A and B 
particles, respectively. The conjugated response fields are denoted by J and c. Stability requires g > 2a f . Green 
functions (correlation and response functions) are obtained by integrating the fields against a weight factor exp(— J). 

The functional J, Eqn. (|l|), possesses the following symmetries under three transformations involving a constant 
continuous parameter a: 



I 

II 
III 



c — > c + a ; 

c-^c + a, t^t + /cv; 
s — » as , s — > a" 1 s , a 



a a 



g^a g 



g -» ag 



Moreover, J is invariant under the inversion: 

IV: c — > — c , c — 
In the particular case a = 0, the time inversion 

V: \fg]gs{x,t) <-> — > 
c(x, t) — » c(x, —t) 



r g/gs{^,- 

c(x, t) - 



-> c(x, — t) — c(x, — t) 



(2) 
(3) 
(4) 



(5) 



(6) 
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yields a further discrete symmetry transformation. The symmetry V distinguishes the KSS from general DP-C 
processes. 

Symmetry I results from the conservation property of the field c. Symmetries III and IV show that dimensionless 
invariant coupling constants and parameters are defined by u — ggp~ e , v — f 2 p~ £ , w = agfp~ £ , and the ratio of the 
kinetic coefficients p = 7/A. Here /i -1 is a convenient mesoscopic length scale and e = 4 — d. 

Dimensional analysis and the scaling symmetry III applied to the Green functions G N ~ ~ = ([sj^fsl^lcj^c^) 
gives 

G N.N;MM = aN ~ NG N,N;M,M^ *}' T ' " ^ a 9> " ^ ^ A ' 7 ' ^ ( 7 ) 

= v N ~ N F N ~. M ~ {{fix, 7/i 2 i}, p~ 2 r, u, v, w, p) (8) 
= (^'"^'^^^^^^O,^ 2 ^,^,,,) , (9) 

where it is assumed that a > for convenience. Of course, Eqn. (||) cannot be used if a = 0. 

The critical scaling properties of the Green functions can be extracted from the invariant functions F and F' 
by applying the renormalization group. To extract UV-finite quantities from the field theory one introduces bare 

o 1/2 

fields and parameters and renormalizes by appropriate Z-factors. For example one uses the scheme s — > s — Z s ' s, 
s — > s — Z^ 2 s, t — > f = Z t t + f c „ f — > / = ZfV 1 l 2 p e l 2 , etc. Here f c denotes the critical value of f. The Z-factors 
have to absorb all the UV- infinities (the e-poles in dimensional regularization) . They can only depend on the invariant 
parameters u, v, w, and p. 

The objects of the calculation are the vertex functions r~ ~ , i.e., the one-particle irreducible amputated 

diagrams with N s-legs, N s-legs, M c-legs, and M c-legs. It is easily seen that diagrams with loops do not contribute 
to vertex functions with M > 1. Thus, these vertex functions are trivial and given by the corresponding terms 
displayed in the dynamic functional J, Eqn. (0). Hence, the renormalizations are trivial: 

c — c, c=c, 7 = 7, as — as . (10) 

Symmetry II in connection with the trivial renormalization of c, Eqn. (pj), shows that / is renormalized with the 
same Z-factor as t: Zj = Z T . It follows the simple relation 

At a fixed point different from and 00, t changes according to this relation by a change of the momentum 
scale p — > pi (holding bare parameters fixed) as r — > r(Z) = rl e / 2 . Thus, one finds from the Eqns. (^,^|) the scaling 
properties of the Green functions at a fixed point with finite values for it* , i>* , 1/;* , and different from and 00 (the 
existence of such fixed points can be demonstrated in the e expansion ||^]) 

G({x, <}, r) = l 5 °G({plx, 1 p 2 l 2 t], p- 2 r/l 2 -^ 2 ) . (12) 



In Eqn. (12) Sq = (N + N + M + M)d/2 + (N-q + Nrj)/2 denotes the scaling exponent of G. Sg combines the normal 
and anomalous dimensions of the fields involved in the Green function G. From Eqn. (|l^) one can gather the exact 
values of the dynamical exponent z = 2 and the correlation length exponent v = 1/(2 — e/2) = 2/d, 

The renormalization of quantities invariant under the transformation defining symmetry III, like F or F', involve 
only the product of the field renormalizations Z = Z~ S Z S . Thus, one has a freedom to define one of this factors. With 
respect to Eqns. @ and (|l^) it is convenient to choose the trivial renormalization a = a together with Z s = 1. Then 
the Green functions G have the same scaling properties under renormalization as the invariant functions F . One could 
also define Z s ^ 1 and renormalize a = Z s a. The renormalization properties of F are not affected by this choice. 
Hence one finds the same critical scaling properties of the correlation and response functions as for the previous one. 
It follows the anomalous dimension of the field s as rj = 0. Then the anomalous dimension rj of the response field s 
is given by the logarithmic derivative of Z with respect to the momentum scale p at the fixed point, rj is the only 
scaling exponent that one has to determine by perturbation theory. 

For the KSS processes a vanishes and one cannot follow the strategy of the last paragraph. However, for the KSS 
processes the time inversion symmetry V can be explored. With respect to this symmetry and Eqn. (j^) it is now 
convenient to choose the ratio g/g — g/g trivially renormalized together with Z s = Zi = Z 1 / 2 . Then the Green 
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functions G have the same critical scaling as the invariant functions F' . The logarithmic derivative of Z at the fixed 
point yields in this case the value of 77 = rj. 

In summary, the DP-C processes offer exact relations for some critical exponents including z = 2 and v = 2/d. In 
the KSS case one has rj = rj as in DP, but with another value. In the more general case with u/0 (and af < 0) one 
finds 77 — 0, which yields via the relation f3 — v(d + rj)/2 the exact order parameter exponent /3 = 1. The discrepancy 
with the exact and the simulational result for v by roughly 10% may have the origin in corrections to scaling. It is 
therefore desirable that the authors reconsider their simulations and provide a careful analysis of such corrections. 
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